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Linear Time-Invariant (LTI) Systems with Random Inputs 


10.2.2 Linear Time-Invariant (LTI) Systems with Random 
Inputs 

Linear Time-Invariant (LTI) Systems: 

A linear time-invariant (LTI) system can be represented by its impulse response (Figure 10.6). 
More specifically, if X(t) is the input signal to the system, the output, 7(0, can be written as 

r oo poo 

7(0 = J _ rJ h{o.)X{t - a) da = J _ x X(a)h(t - a) da. 

The above integral is called the convolution of h and A, and we write 

7(0 = h(t) * X(t) = X(t) * /?(/). 

Note that as the name suggests, the impulse response can be obtained if the input to the system 
is chosen to be the unit impulse function (delta function) x{t) = S(t). For discrete-time systems, 
the output can be written as (Figure 10.6) 

Y(n) = h(n ) * X(n) = X(n) * h(n) 

OO OO 

= X h(k)X(n ~k)=Yj X(k)h(n - k). 

k= —co k= -co 

The discrete-time unit impulse function is defined as 

.1 n = 0 

§{n) = | 

10 otherwise 

For the rest of this chapter, we mainly focus on continuous-time signals. 
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X(t) 


* h(t ) 


LTI System 


- *Y(t) 

Y(t) = h(t) * X(t) 


X(n) 


■> h(n ) 


LTI System 


- *Y(n) 

Y(n) = h(n ) * X(n) 


Figure 10.6 - LTI systems. 


LTI Systems with Random Inputs: 


Consider an LTI system with impulse response h(t). LetX(0 be a WSS random process. IfX(0 
is the input of the system, then the output, 7(0, is also a random process. More specifically, we 
can write 


7(0 = h(t) * X(t) 

roo 

= I -oo h{a)X(t - a) da. 

Here, our goal is to show that X(t) and 7(0 are jointly WSS processes. Let's first start by 
calculating the mean function of 7(0, //y(0- We have 


MyiO = ^[T(0] 


E 


roo 

I - x h(a)X{t ~ a) da 


r oo 

= J -oo h(a)E[X(t - a)] da 

roo 

= J _ o0 h(a)fi x da 


roo 


= Mx] -oo h ^ da - 


We note that // y(0 is not a function of t, so we can write 


roo 

My(t) = ji Y = Hx) -csM a ) da. 

Let's next find the cross-correlation function, R XY it\, t 2 )- We have 
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i>oo 

Rxfti’ h ) = E \. x i.h) Y ^2)\ = E ml-JWXfi - a ) da 

" poo 

= E J ryD h(a)X(t ] )X(t 2 ~ a) da 

poo 

= J _ 00 A(a)£[X(? 1 )LSr(?2 “ a )l da 

TOO 

= I _ oo h(a)R^{t j, “ a) 

poo 

= J _ co h(a)Rj i (t l ~ t 2 + a) da (since X(0 is WSS). 

We note that R^t^ t 2 ) is only a function of r = /, - / 2 , so we may write 

roo 

e Xy( t ^ = J -oo h(d)Ry(x + a) da 

= Kx) * R£- x) = h{~ x) * ^t). 

Similarly, you can show that 

Ryir) = h(r) * h{ - r) * R^x). 

This has been shown in the Solved Problems section. From the above results we conclude that 
X(t) and 7(0 are jointly WSS. The following theorem summarizes the results. 


Theorem 10.2 

LetX(0 be a WSS random process and 7(0 be given by 

7(0 = h(t) * X(t\ 

where h(t) is the impulse response of the system. Then X(t) and 7(0 are jointly WSS. 
Moreover, 

1. n y(0 = ju Y = ju^ x . o0 h(a) da; 

2. Rxy(x) = //( - 0 * 0 = l°^ 00 //( - a)R x {t - a) da; 

3. Ry {0 = h(r) * h{ - r) * R^(r). 

Frequency Domain Analysis: 

Let's now rewrite the statement of Theorem 10.2 in the frequency domain. Let H(f) be the 
Fourier transform of h(t), 
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m = nm = \Zj(t)e~ 2 ^ dt. 


H(f) is called the transfer function of the system. We can rewrite 

roo 

Ry~ V'X]-JW da 


as 


R y 

Since h(t) is assumed to be a real signal, we have 

F{h{-t)}=H{-j)=H*(j), 

where * shows the complex conjugate. By taking the Fourier transform from both sides of 
Rx&) = R%(j) * h( - r), we conclude 

S X y(f) = S^/)/I( ~f) = S^j)H*(j). 

Finally, by taking the Fourier transform from both sides of Ry(z) = h(z) * h( - z) * R^r), we 
conclude 


Syif) = S]£f)H* (f)H(f) 
= S x <f)\H(J)\ 2 . 

S,0 = Syj)\H{f)\ 2 


Example 10.14 

Let X(t ) be a zero-mean WSS process with R^z) = e ~ I T L X{t) is input to an LTI system with 


i m 


y/l+47r 2 / 2 

0 


I/I < 2 


otherwise 


Let 7(0 be the output. 

a. Find Ry(t) = E[Y(t)\. 

b. Find Ry(z). 

c. Find E[Y(t) 2 ]. 
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• Solution 

o Note that since X(t) is WSS, X(t) and Y(t) are jointly WSS, and therefore Y(t) is 
WSS. 

a. To find u y(t), we can write 


My ~ Mxm 

= 0-1 = 0 . 

b. To find Ryir), we first find .S’//). 

Syif) = S^J) | H(f) | 2 . 

From Fourier transform tables, we can see that 

2 

" 1 + (2tt/) 2 ‘ 

Then, we can find Sy(f) as 

Sy(f) = S^f)\H{f) | 2 

( 2 I/I <2 

10 otherwise 

We can now find R /r) by taking the inverse Fourier transform of Sy(f). 


where 


c. We have 


Ryir) = 8sinc(4r), 


sinc(/) 


sin//) 


= Ry(0) = 8. 


previous 
next —» 
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